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1 INTRODUCTION 



Branching stochastic processes have always been very interesting for mathematicians and 
physicists. They describe well a multitude of phenomena from chain reactions to populational 
dynamics. On the other hand, ordinary |^ and multi-scale ^ percolation play a crucial role 
in many applications 0]. In the work (see also 0) on the basis of the idea of recoding |^ 
necessary and sufficient conditions for the discrete hierarchical (multi-scale) model of fracture 
(percolation of defects) were obtained. 

In this paper we give an account of matter and (basically) a solution of a new class of 
problems synthesizing percolation theory and branching diffusion processes. Such problems 
arose naturally in the recent investigations of the global geometry of inflationary early Universe 
(see and references therein). 

In Section ^ we present a general formulation of problem of inflnite-scale percolation in 
a new type of branching diffusion processes. The novelty is the diffusion not in the real space 
where branching "particles" live, but in the space of parameters distinguishing the "particles" 
each other. 

We do not specify in the following, whether the A-diffusion on a (i-dimensional manifold, 
which we assume from the very beginning, corresponds to some stochastic differential equation 
(SDE) in the sense of Ito or in the sense of Stratonovich. In the applications P| we use the 
Stratonovich's symmetrized calculus because white noise there is the limit of a colored one 
with small time correlation. The generating operator (A-operator) of the diffusion in local 
coordinate frame looks like ^ {i,j,k = 1,2, ... ,d ; the summation over repeated indexes is 
assumed) 

i) after Ito 

AifiX) = la^\X)a'^{X)^j^f{X) + r{X)^f{X) (1) 

ii) after Stratonovich 

Asfix) = IcT'Hx)^ + ^W^fiX) (2) 

and corresponds to SDE (in the proper sense) 

dXi = r{X) dt + a^'^(X) o dW^ (3) 

Here XI G R^, is the d-dimensional Wiener process. We add to the random walk (^ 
a branching with intensity n{X) (i.e. the probability density of branching of the "particle" 
at X during the time At is n{X)/S.t ). Then we relate with such a branching diffusion some 
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process of breaking up of D-dimensional cubes and their coloring and obtain an infinite-scale 
(in the t oo limit) percolation problem. 

After some preliminary (Section |^) the solution of the problem is stated in Section ^ where 
necessary and sufficient conditions for percolation are obtained. In Section |^ we illustrate our 
method in the case of the simplest model. In the Summary (Section ^ possible generalizations 
are discussed. 



2 THE PROBLEM 

Let the basic (i-dimensional diffusion process generated by the operator A (|l]), (|^) given 

x, = x] = Af{X) (4) 



lim E 



\ At 



where G A^, is a ci-dimensional manifold. An action of A is endowed by the boundary 
condition 



fiX)\xer = (5) 

where T G A4 is a closed subset of A4 (an absorbing boundary). 

On A4\T smooth functions n(X), JF*(X), o"*''(X) are defined (see (^), where n{X) and 
a^^{X) are assumed to be positive, such that F is accessible. Let in each branching instead 
of one "particle" at point X appear the given number r of "particles" at that point, which 
continue to evolve as branching diffusion process independently from each other. When a 
"particle" reaches the absorbing boundary F, it stays there forever without branching. 

Now we will associate this branching diffusion with the following picture. Let at t = a 
Z)- dimensional cubic net be given, consisting of cubes of unit size. With each cube of this net 
independently we set in correspondence a random point Xo G with the probability density 
^{Xq). The subsequent evolution of each cube is independent from the rest of the net. This 
evolution is determined by the abovedescribed basic branching diffusion on A^. At t = r 
trajectories X^^'^ (z = 1, . . . , r) of random walkers (^ start at Xq. We divide the unit cube (call 
it cube of zero level) into r smaller equal cubes of first level (we assume that r = k^, so that 
the cubes of first level are k times smaller than the unit cube). Each first level cube is set in 
correspondence with one of the points xj:^^ — the points where the z-th "particle" branches for 
the first time. If the "particle" was absorbed at F, we attribute to the corresponding first level 
cube the proper point Xi E T, color it black and leave it in peace. We divide the other cubes, 
which should not be colored this time (we call them "white" or "living" cubes), once more 
into the r = k^ cubes of second level and do the same procedure, starting from X^^ G A1\F. 
Thus, coloring some cubes at each level q (when Xjj^^ G F ) and dividing white cubes further we 
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will obtain some infinite-scale (provided that the process does not degenerate at finite level) 
picture of black cubes of different sizes in a "sea" of living white cubes (see Fig. |I|). 



Denote Bt the multitude of all black cubes of all sizes (and Wt — of all white cubes) 
belonging to a particular unit cube at the time t. We assume that two cubes of (possibly) 
different sizes are connected if they have common face. Now we are ready to formulate the 

Problem: Consider the whole net of (zero level) cubes and denote by the union of all the 
multitudes Bt belonging to all unit cubes in the original cubic net. Let j{X), n{X) and 
A be given. Does Bt percolate in the t — >■ cx) limit (we will denote the corresponding 
multitudes of cubes as Boo and Boo respectively)? 

Note: Percolation of Boo means that a connected non-selfintersecting path a. s. exists on Boo 
from some point (call it origin of the coordinate frame) up to infinity. 



First of all, let us elucidate the condition of non-degeneracy of the process of breaking up of 
the cubes. Let fi{t, V) be the number of random walkers within the region V G A4\T at the 
time t. Introduce the generating function 



Noting that the evolution of u{t, X, z) is driven by the diffusion of particles and by the branch- 
ing processes, we obtain 



u{t + /\t,X,z) = {l-n{X)At)E{u{t,XAt,z)\ Xq = X} + n{X) At {u{t, X, z)Y (7) 



Here we considered the whole tree of the branching diffusion process from t = to t + At 
and divided it into two parts — a part from t = to t = At and a part from t = At to t + At. 
Up to the first order in At, we have only two mutually exclusive possibilities of evolution of 
the process from t = to t = At. 

1. The single original particle does not branch during this time and it diffuses to new place 
X/^f The tree of the branching diffusion process from t = At to t + At differs from the 
one corresponding to the interval of time [0, t] only by its origin X/^f This possibility is 
represented by the first term in the r.h.s. of Eq.(|^). 

2. The original particle branches just one time during the period from t = to t = At. 
There appear r species of the tree of the branching diffusion process and, in the first 
order in At, we should not distinguish them from the tree corresponding to the interval 
[0,t]. This possibility is represented by the second term in the r.h.s. of Eq.(0). 



3 PRELIMINARY 




(6) 
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Then, using Eq.(|D, one easily derives from (|^) the (backward) differential equation for the 
generating function 

X,z)=A u{t, X, z) + n{X) {u\t, X, z) - u{t, X, z)) (8) 
and the boundary and initial conditions 

u{t,X,z)\x^r = l (9) 

u{0,X,z) I otherwise ^^^^ 

As a consequence of the definition (P) we can write down the equations for the (factorial) 
moments of fi{t, V) differentiating the Eqs. (p|)-(p!0|). Note that 

mh{t, X) = E{fi{fi-l)...{fi-l)\Xo = X} = ^u{t,X, z)l=, (11) 

and 

mi{t,X) = E{fi\t,V)\ Xo = x} = [z§^)'u{t,X,z)U=i (12) 

In particular, the equation for the first moment mi{t,X) = z^u{t, X, z)\z=i (the average 
number of random walkers in V at the time t ) is 

fmi(t,X) = Ami{t,X) + {r - 1) n{X) mi{t, X) (13) 
m,{t, X)\xer = (14) 

Eqs. (|l^)-(0) have an asymptotic {t oo) solution 

mi{t,X) (xe^^'MX) I dYniiY) (16) 
Here ipi{X) is the unique strictly positive real eigenfunction of the marginal problem 
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i^i(X) + (r - 1) n(X) ^i(X) = AiVi(X) 



(17) 



^i(X)Uer = (18) 

and Ai is the corresponding (real) eigenvalue. The function 7ri(X) (invariant density) is the 
unique strictly positive eigenfunction of the adjoint equation with the same eigenvalue Ai 

it 7ri(X) + (r - 1) n{X) 7ii{X) = Ai7ri(X) (19) 
7ri(X)Uer = (20) 

The normalizations follows 

/ Tri{X)dX = l; [ tti{X) ^i{X) dX = 1 (21) 
JM\r JM\r 

One can see from (|16D that if Ai < 0, then the branching process a. s. degenerates in the 
limit f — i> oo (in our notations Woo = 0)- If Ai > 0, then the branching process is supercritical 
and Woo 7^ • We assume the latter case in the rest of the paper. 

Let us introduce two other useful functions. The first one is the probability of the event 
that the particle gets absorbed at T at some time less or equal t provided that it started to 
diffuse from some point X G A^\r at the time t = and didn't branch before being absorbed. 

p(X,t) = P{Kt|Xo = X} (22) 

Consider the moment of time t + At. It is obvious from the definition above that (in the first 
order in At) 

p(X, t + At) = {l- n{X) At) E { p(XAt, t) | Xq = X} (23) 

and the following (backward) differential equation as well as the initial and boundary conditions 
follow immediately from ( pSj) and (^) 

|p(X, t) = A p{X, t) - n{X) p(X, t) (24) 
p(X,t)Uer = l (25) 
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In the rest of this paper we will use only the stationary probability p{X) of absorption 
of a particle, which starts its random walk at X, without branching before absorption. It is 



given by the stationary solution of Eqs. (p^)- (l26|) . Since n(X) is positive and the maximal 
eigenvalue of A with zero boundary condition on F is negative, A = is not an eigenvalue of 
{A — n{X)) and the stationary equation 

Ap{X)-n{X)p{X) = ■ p(X)Uer = l (27) 

has unique solution. Of course, the relation limj^oo p{X, t) = p{X) is satisfied. 
Then, let us introduce the probability density in Y 

K{X,Y,t) = P{^,{Y)\Xo = X} (28) 

where ^t{Y) is the event that the particle branches for the first time in the infinitesimal volume 
Y + dY at some time less or equal t, provided that it starts to diffuse at X at t = 0. Consider 
the evolution of this quantity after time At in the first order in At. One has 



K{X,Y,t) = n{X) M 5{Y - X) + {I - n{X) /\t)£{K{X^uY,t + ^t)\ Xq = X} (29) 

Here the first term represents the event that the particle branches during the time interval 
[0, At] and we should not distinguish the position of the particle at t = and t = At in the first 
order in At. The second term represents the event that the particle does not branch during 
the time interval [0, At] and in the first order in At the only difference we should take into 
account for that particle's trajectory during the time interval [At, t + At] is its random initial 
point at the beginning of that interval. These are the only two mutually exclusive events 
which exist in the first order in At. From the definition ( pSj ) and Eq. (^) one easily obtains 
the (backward) differential equation and the boundary and initial conditions 

|if (X, y, t) = AxK{X, y, t) - n{X) K{X, r, t) + n(y) Six - Y) (30) 
K{X, Y, t) \xer = K{X, Y, t) |yer = (31) 
K{X, Y, 0) = 6{X - Y) (32) 
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Here Ax denotes the generating differential operator of diffusion acting on the first 
variable of K{X,Y,t). In what follows we will need only the stationary probability density 
K{X,Y), which satisfies the equations 

AxK{X, Y) - n{X) K{X, Y) + n{Y) 6{X -Y) = (33) 

K{X,Y)\xeT = K{X,Y)\YeT = (34) 

One can write down the expression for K{X, Y) through the complete orthonormal set 
of eigenfunctions Ks{X) and eigenvalues Ks of the marginal problem 

AKs{X) - n{X) K,{X) = k,K,{X) (35) 

Ks{X)\xer = (36) 
Recalling the following properties of the eigenfunctions 



/ K,{X)K,{X)dX = 6,j ■ Y.Ks{X)Ks{Y) = 6{X-Y) (37) 
one easily obtains 

K{X, Y) = -n{Y) J2 -Ks{X) K,{Y) (38) 

s=l '^s 

Now, with the functions p{X) and K{X, Y) in hand, we are ready to solve the problem 
stated in the preceding section. 



4 INFINITE-SCALE PERCOLATION 

In this section we will find the percolation characteristics of Bqo using the renormalization 
relations, which arise as a result of the recoding procedure. This procedure was introduced in 
the case of the simpler model of discrete hierarchical fracture in [0] and was investigated in 
details in |^ . 

Consider the multitude Bt , belonging to a particular unit cube, at some large time t. 
Let us concentrate for the moment on the percolation characteristics of this unit cube. Let us 
define the notion of "strongly defective" cubesQ. 

■^We call black cubes "defects" to resemble the rock fracture theory m. 
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Definition 1 A cube of level L is called "strongly defective" of rank (SD^^ ) if it is black. It 
is called "strongly defective" of rank m > (SD^"*), if either it is black or the following two 
conditions are satisfied: 

1. There exists a connected cluster o/SD^]^^"* cubes of level L + 1 which connects each pair 
of faces of the given cube; 

2. Strictly more than a half of the surface of each face of the cube is covered by SD^^^^ 
cubes belonging to one of the clusters defined in 1) 

A cube is said to be SD if it is SD*-'^-' for some n > 0. 

We assume in the main body of the text A; > 3. In the k = 2 case more care is needed 
due to the fact that there are no nontrivial configurations satisfying condition 2 (see Section 
P). The meaning of this definition is revealed by the following obvious 

Proposition 1 If any number of SD cubes of some level L are attached face to face, then there 
is a connected cluster of black cubes (consisting, maybe, of cubes of different levels K > L) 
running through them. If a cube is SD^^^^ then it is SD*-"^ for every n> m. 

Note, that it is the condition 2 of the Definition |l| that ensures the percolation of the 
black color through the attached faces of SD cubes. Thus, SD cubes are as good as black ones 
when we deal with the percolation of the black color through the net of cubes. It is this kind 
of substitution of black cubes by SD cubes that was called "recoding procedure" in [0 and 
that enables us to follow the "renormalization approach" of 0, 0]. 

An opposite proposition (that SD cubes always consist of SD cubes) is false (see Fig.|]). 
That's why we can obtain only sufficient (but not necessary) condition for percolation assuming 
that the probability of level L cube to be SD can be computed neglecting such configurations 
as in Fig.|^ 

P. {SD J . P. {B., . (1 - P. (B.,) P. { lef ) (^«) 

where the subscripts L denote the level of the cubes and the subscripts X take into account 
the fact that we deal with probabilities depending on the value of the parameter associated 
with the given cube. We are interested in calculation of Px {SDq}. Let us remind that SD 
means SD(") for some n > 0. In other words, we should find the limit 

P{X} = JimP^{SDS'^)} (40) 
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Let us forget for a while about the time dependence of the problem and use the stationary 
equation 

Px {SD^")} = p{X) + (1 - p{X)) j^^^ KiX, Y) Qo,k{PY {SD^"-^)}) dY (41) 

where p{X) and K{X,Y) were introduced in the preceding section (see (^), (^)-(^)). 
Equation (^Tj) gives the probability for the zero level cube to be SDq"'' in terms of the prob- 
abilities of the first level cubes to be SDj^""^'' provided that the parameters associated with 
them are such that the "particle" corresponding to the zero level cube is the "parent" of the 
particles corresponding to the first level cubes. 

The combinatorial function QD,k{p) which appeared in (^Tj) is a polynomial in p of r-th 
degree. It counts the probability of SD configurations of level L + 1 cubes inside of cubes of 
level L (we assume for the moment that p is the probability of black color and q = 1 — p is 
the probability of white one). We will consider in this paper only the simplest case of two 
dimensional cubes (squares) and k = 3, although it is only a straightforward combinatorial 
problem to find QD,k{p) for the other casesQ. 

Q2,3(p) =p^ + 9p^q + 20/g2 = 12/ - 31/ + 20p^ (42) 
Equation (^) can be continued by the hierarchy 



Px {SDi"-^)} = p{X) + (1 - p{X)) 1^^^ dYKiX,Y) QnAPy {SBf^}) (43) 

Px {SD^i} = p{X) + (1 - p{X)) J^^^ dYK{X,Y) QnAPy {SDf }) (44) 
Then, recalling that 

Px{SDS,,e.e/} = PW (45) 

we easily find 



>x {SD^")} = F[F[. ..F[p{X)] ...]] = F4p{X)] (46) 



n times 



"^One should remember, however, that computational difficulties grow exponentially fast with growing D or 
k, and the finding of the functions QD,kip), although remaining a finite combinatorial problem, becomes very 
cumbersome. 
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Here denotes the n-th iteration of the integral operator 

F[ct>{X)] = p{X) + (1 - p{X)) j^^^ dYKiX, Y) QnA<t>iY)) (47) 

From Eqs. pT|)-(^7D one can derive that in the n — oo (t — ^ oo) hmit the iterations 
converge (under some restrictions on p{X) and K{X, Y) which are necessary to guarantee 
the uniqueness of the fixed point) to the solution of the following integral equation 

P {X} = F[P {X}] (48) 



It is worth noting that after Eqs. (|41| ) - (|4^) it becomes clear why we can substitute a t- 
dependent percolation problem by n-dependent hierarchical one. Indeed, up to now we didn't 
mention that not all "white" cubes of given level break up simultaneously and consequently 
there exist many "white" cubes of different sizes at the time t. However, our approach is 
insensitive to this difficulty because different cubes can join the "renormalization flow" at 
different scales, but all of them approach eventually the same fixed point (0). The only 
modification that can arise from this notion is the somewhat stronger condition on the limiting 
rate t — > oo when one proves the coincidence of lim^^oo Px {Bt} and P {X} from (40). 



In some cases, when the full transition probability density K(X, Y) ( pSj) can be approx- 
imated by a finite sum like 



niY)j:-K,iX)K,{Y) (49) 



the corresponding nonlinear integral equation ( |48| ) reduces to an algebraic one. This fact can 
be useful in applications [§. 

Suppose, that we found the solution of (^). Then the solution of the percolation problem 
is straightforward. The quantity 



p,= / 7(Xo)P{Xo}dXo (50) 
JM\r 

has the meaning of the probability that the cube of zero level is SD. If p^'' is the percolation 
threshold for the site percolation problem on a cubic Z)-dimensional lattice Z^, then, relating 
black sites to SD cubes, white sites to non-SD cubes and recalling the Proposition |^, we come 
immediately to 



Theorem 1 The sufficient condition for percolation of Boo is 

P* > pff (51) 
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The necessary condition for percolation is the negation of the sufficient condition for non- 
percolation. Let us introduce the notion of "strictly closing" (SC) configurations to obtain the 
sufficient condition for non-percolation. 



Definition 2 A cube of level L is called "strictly closing" of zero rank (SC]^ ) if it is white. 
A cube is called "strictly closing" of rank m (SCj^'*) if it is either white or it consists of such 
configuration ofSCi^i^^ cubes that there is no connected cluster of the remaining cubes of level 
L + 1 belonging to SC^-* cube which could connect any pair of faces of that cube. A cube is 
said to be SC if it is SC^"-* for some n > 0. 

One can easily verify the 

Proposition 2 Any number of SC cubes, attached face to face, do not contain any percolating 
path of black cubes. If the cube is SC*-™-*, then it is SC^"-* for every n> m. 

Now it is clear that our "renormalization" procedure can be applied. Let us denote by 
the bar all quantities concerning the SC cubes. In complete analogy with the Eqs. (p9|)-(|48|) 
we obtain 

Px{sc(")} = (l-p(X)) j^^^ dYK{X,Y)Qj,^,{Px{sct''^]) (52) 
Px {SCi^li} = $[£[^[Px {SCf }] ...]] = - p{X)] (53) 

n times 

where we used Px {^^^'^Ulyievei] — ^ ~ P(^) ^^^1 the definitions 

<l>[0(X)] = (1 - p{X)) j^^^ dYK{X, Y) QnMY)) (54) 

P{X} = JimPx{SC('^)} (55) 
Then, we have (under the condition that $[0(X)] has a unique fixed point): 

P{X} = <|.[P{X}] (56) 
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Here Qu kil) is an analog of QD,k{p) for closing configurations {q = 1 — p is the probability 
of white color, p is the probability of black one). We will use in the next Section only the 
simplest combinatorial function of that type 

^2,3 (9) = q^ + 5q^P + lOg^p^ + 4gV + = 3g^ - 7q^ + W + q^ (57) 

P {X} is the probability for the zero level cube to be SC. Associating white sites of black 
site percolation problem on to SC cubes, black sites to non-SC cubes and keeping in mind 
the Proposition |^, we come to 

(2) (2) 

Theorem 2 Let p'^' = 1- q'^' he the non-percolation threshold for the black site percolation 
problem on the . Denote 

q,= f 7(Xo) P{Xo}rfXo (58) 

JM\T 

Then the necessary condition for percolation (negation of the sufficient condition for non- 
percolation ) of Boo is 

q* < q^S (59) 



5 THE SIMPLEST MODEL 

Here we will investigate the simplest model of the type described above. It corresponds to 
homogeneous diffusion on the line segment M. = [0;/];r = {0;l};d = 1 with 

jF(x) = ; cr(x) = a = const (60) 



n{x) = 0? = const ; 7(0;) = - (61) 

t 



^/(^) = Y^/(^) (62) 



and to D = 2; k = 3;r = 9 case of the cubic net. From Eqs. ([T7|) , ([iq) we obtain 



2^2 



Ai = - ^ (63) 
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and the non-degeneracy condition is ^ > |. Assuming this, one can derive from equations 



cosh(iC — z) 
cosh(iC) 



(64) 



K{z,u) 



y-£^sinh(C-2;) sinh(M), ii z > u 
y^j^sinh(C-M) sinh(2;), if z < u 



(65) 



where we have introduced the dimensionless variables and the parameter (, which is in fact the 
only dimensionless parameter in this model and on which the behaviour of the model depends 
essentially 



\f2ax \/^ay \plal 

z = ; u = ; C = 

a a a 



(66) 



Equation (H3) now reads 



P{z} 



+ sinh(2;) / 

J z 



cosh(|C — z) I cosh(|C "~ -2) 

cosh(|C) ^\ cosh(iC) 

C sinh(C- w) (o i \\ 
du — r^^-^ — Q2,3(P 



sinh(C) 



sinh(C-^) r du^^^O,,,{P {u]) 
Jo smh[Q 



(67) 



Introduce the function 



, sinh(tt) ^ r N 
smh(C) 



(68) 



Because we have p{( — z) = p{z), symmetry reasons lead us to P {C — -z} = P {z} and Eq. (|67|) 
now reads 



One can derive an ordinary differential equation for P{z) from (|6^) , (|69D , but it is too 
involved and we don't produce it here. Instead, we present the result of the computer solution 
of inequality (|5T|) in terms of which was calculated using (|68D , (|69|) and assuming the known 
value of p^2^ = 0.59. Thus, the sufficient condition is 
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C < Csuff = 3.29 



(70) 



Equations for the necessary condition are (see 



P{z} 



1 - 



cosh(|C 



cosh(|C) 




sinh(C — z) g{z) + sinh(z) g{( — z)) 



(71) 







sinh(C) 



(72) 



Our estimate for the necessary condition ( |55| ) is too crude for this particular model and 
we don't produce the result here. However, we think that the sufficient condition ( ffOD is 
reasonably close to the exact critical value of (. This can be confirmed by calculation of the 



next level corrections to it (compare with Eq. (|73D, (0), where such corrections are the leading 
ones) . 



6 SUMMARY AND DISCUSSION 

We have described a new class of problems, synthesizing ones from familiar percolation and 
branching diffusion. To find the percolation characteristics, one has to solve the analytic 
relations (13)- (H) and (il)-(|5|). 

Let us emphasize some possible generalizations. First of all, it is worth noting that face- 
connectedness, which we assumed throughout this paper, can be replaced by a great variety of 
other "quasilocal" definitions of connectedness. By this we mean a definition, which would deal 
with some finite clusters of black cubes (possibly of different sizes) with some rules restricting 
their configurations. If something like Propositions |l|, ^ holds, then only minor modification of 
our approach, reducing mainly to change of combinatorial functions QD,k{p) ^"^^ QD,ki.l) ^^^d, 
maybe, the number of levels of hierarchy involved in construction of recurrent formulas of type 
(p!D, (|5^) should be expected. For example, one can consider two black cubes to be connected 
iff the minimum path from one to the other contains no more than Rq white cubes. This would 
correspond, in some sense, to defects with "spheres of influence", which were investigated in 
some one-scale models (i?o = corresponds to the face to face connectedness, i.e. the 
present paper). With this modiflcation, the problem will have an additional parameter Rq 
and, presumably, non-trivial dependence on that parameter. One can consider a new problem 
taking the other parameters of the model flxed and letting Rq vary. It is clear that even if 
there were no percolation of this "spheres of influence" for small i?o, it would presumably occur 
for safficiently large "radius of sphere" i?o, and, therefore, we would have a percolation phase 
transition at some intermediate value of Rq. 

However, the consideration of more than two levels of hierarchy appears to be necessary 
already for i?o = if we consider k = 2 case. Here we present the analogs of the equations (P7| ) 
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and ( ^41) for the k = 2 case. Because at the first step of breaking up there are no nontrivial SD 
(and SC) configurations, one has to consider two steps of breaking up (and, therefore, three 
levels of hierarchy) in order to derive the corresponding integral operator. The form of this 
integral operator is given by 



F(2)[P {X}] = p{X) + (1 - p{X)) J dYK{X, Y)[p\Y) 
+ p3(r)(l - piY)) J dZKiY, Z)Qi'i{P {Z}) 
+ p2(r)(i - p(y))2 J dZdUK{Y, Z)K{Y, [/)Qg(P {Z}, P {U}) 
+ p(y)(l - p{Y)f J dZdUdVKiY, Z)K{Y, U)K{Y, V)Q%{P {Z}, P {[/}, P {V}) 
+ (1 - p{Y)Y J dZdUdVdWK{Y, Z)K{Y, U)K{Y, V)K{Y, W) 

xQ^^^iP{Z},P{U},P{V},P{W})] (73) 

where the combinatorial functions Q2^2{Pi)j Q'^Mp^^P'^)^ ^2^2(^1 1^2,^3), ,P2,P3,P4) are 

introduced which count the "strongly defective" configurations of black cubes of two subsequent 
levels L + 1 and L + 2 of hierarchy within the "parent" cube of level L for the cases when 
only one (two, three or four, respectively) of the cubes of level L + 1 is broken up to cubes of 
level L + 2 (and the first term in the integrand takes into account the only SD configuration 
when none of L + 1 level cubes is broken up). Their expressions as polynomials of pi are too 
involved and do not seem to contain any clarifying information, so we do not produce them 
here. Equations (148|)-(|5TD are implemented with instead of F[-] . Analogously, for the 



necessary condition we should use in ( |56D the following expression instead of $[■] (see (Q)) 

$(2)[p {X}] = (1 - p(x)) I dYKix, y)[(i - piY)y 

X J dZdUdVdWK{Y,Z)K{Y,U)K{Y,V)K{Y,W) 
xQ'-^1{P{Z},P{U},P{V},P{W})] (74) 

with obvious notations. 

Second, one can generalize the branching diffusion part of the problem. For instance, 
the number r can be considered to be random (possibly as a function of ). This, of course, 
will make the problem of definition of "strongly defective" and "strictly closing" configurations 
more complicated because in this case we should deal with face to face attached cubes which are 
broken up in substantially different ways and we should pay much more attention to ensuring 
the existence of something like Propositions 0, ^. However, we think that the revision might 
be only technical (taking into account that the branching diffusion part of the problem may 
be generalized to this case easily [|10| ) and in physical applications such a modification can 



always be absorbed by some corrections to the function n{X). 
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One can improve the evaluation of both necessary and sufficient conditions considering 
three (or more) steps instead of two. In principle this can be done and one would obtain the 
equations similar to (|73D and (0), but computational difficulties increase very fast making 
this method hardly applicable. 

For applications it is very interesting to develop any method of approximate solution 
of equations (|48|) , (|56|) . Sometimes, approximations, which reduce K{X,Y) to degenerate 
kernel, are not so bad. 

Anyway, the investigated problem has physical significance [§] . It may occur that slightly 
complicated models of fracture, currently used in geophysics P -[|^ as well as some general- 
izations of the models currently used to describe the intermittent behaviour of high-energy 
scattering processes also will lead to similar problems. We hope that the applications will 
not be limited to these known areas. 



A.M. is indebted to A.D.Linde for suggesting the physical problem [Q, that led us finally 
to the present work. He is also grateful to A.A.Starobinsky, Ya.G. Sinai, L.Bogachev and 
O.Filippova for useful discussions on various aspects of this work. We highly appreciate the 
careful and kind reference and are grateful to the referee for very useful suggestions concerning 
especially the definitions of SD and SC cubes. 
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